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Abstract: We investigate the realization of two bouncing paradigms, namely of the super¬ 
bounce and the loop quantum cosmological ekpyrosis, in the framework of various modified 
gravities. In particular, we focus on the F{R), F{G) and F{T) gravities, and we recon¬ 
struct their specific subclasses which lead to such universe evolutions. These subclasses 
constitute from power laws, polynomials, or hyper geometric ansatzes, which can be ap¬ 
proximated by power laws. The qualitative similarity of the different effective gravities 
which realize the above two bouncing cosmologies, indicates that a universality might be 
lying behind the bounce. Finally, performing a linear perturbation analysis, we show that 
the obtained solutions are conditionally or fully stable. 
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1 Introduction 


According to many theoretical arguments and observational indications, the universe exhib¬ 
ited an early accelerating phase called inflation [1, 2], The inflationary paradigm, amongst 
others, is very efficient in solving the flatness, horizon and monopole problem, and offers a 
consistent mechanism for the generation of the primordial fluctuations or primordial gravi¬ 
tational waves [3, 4] (see also [5, 6]). Hence, current observational research [7, 8] is focused 
on revealing the Universe’s evolution during all eras and theoretical research is focused on 
embedding all the striking new results in a consistent theoretical framework, and if possible 
under the same theoretical framework. 

One alternative to the standard inflation description of the early acceleration is pro¬ 
vided by bouncing cosmologies [9, 10], in which one avoids the appearance of the initial 
singularity. Some of the most appealing representative theories in this paradigm, are the 
ones that use scalar fields [9-22], the ones that use various models of modified gravity 
[23-35], and the ones that use the Loop Quantum Cosmology [36-44] (LQC) theoretical 
framework. In most of these scenarios, the bounce realization relies on the existence of 
matter helds, with an equation of state of suitable form in order to achieve the necessary 
violation of the null energy condition [9-12]. Bouncing cosmologies may have primordial 
instabilities, which make the contracting phase a rather problematic era [13, 14], however 
one can solve this problem in many ways, as for instance in the framework of ekpyrotic 
scenario [11, 12]. For relevant works on bouncing cosmologies see [45-47] and for a study on 
the non-existence of open or flat bouncing cosmologies in Einstein gravity, see [48]. Finally, 
note that there might be observable signatures of the bouncing phase [49-51]. 

In the present work we are interested in investigating the superbounce and the LQC 
ekpyrotic scenarios in the framework of modified gravity [52-54]. Since modified gravity 
is one of the two main directions in which one can achieve also the late-time acceleration 
(with the other one being dark energy), the above incorporation of bouncing solutions 
in modified gravity becomes more important under a unified picture, namely to be able 
to obtain the bouncing-initial phase of the universe as well as its late-time acceleration, 
simultaneously. In particular, we will investigate the above bouncing behaviors in the 
context of F{R) [55-71] (for recent reviews see [72, 73]), F{G) [74-83, 86, 87] and F{T) 
[88-113] theories. The techniques that we shall use in order to provide a modified gravity 
description of the bouncing cosmologies are quite renowned [65, 66, 114] (see also [67-69]). 
For related works on bounce cosmology reconstruction from F{R) theories see [115, 116], 
and for F{R) ekpyrotic cosmology see [117]. Finally, for completeness, we investigate 
whether the obtained solutions are stable under linear perturbations. 

The manuscript is organized as follows: In section 2, using very well known reconstruc¬ 
tion techniques, we investigate the F{R) modified gravity, and in particular we extract the 
F{R) forms that give rise superbounce and LQC ekpyrotic evolution. For the superbounce 
case we provide a complete analytic solution in closed form, while for the case of LQC 
ekpyrosis we extract the large and small curvature limits. In Section 3 we perform the 
same analysis for the case of F{G) gravity, and similarly, in Section 4 for the F{T) grav¬ 
itational paradigm. In Section 5 we investigate the stability of the obtained F{R), F{G) 
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and F{T) solutions. With regards to the F{R) and F{G), the stability analysis is focused 
on examining the dynamical system formed by the resulting FRW equations, under linear 
perturbations. Additionally, the resulting conditions that ensure stability are thoroughly 
investigated. Finally, Section 6 is devoted to the conclusions. 

2 Superbounce and loop quantum ekpyrotic cosmology from F{R) grav¬ 
ity 

In this Section we will reconstruct the suitable F{R) form that can lead to the superbounce 
and loop quantum ekpyrotic cosmology realizations. Before investigating the reconstruction 
procedure, we briefly review in the following subsection the basic features of F{R) gravity 
in the metric formalism and we describe the geometric background we shall use in its 
cosmological application. For a detailed account on these issues, see [72, 73] and references 
therein. 

2.1 F{R) gravity and cosmology 

Let us review briefly the F{R) gravitational modification and its cosmological application. 
Throughout the article we assume that the spacetime manifold is pseudo-Riemannian, 
locally described by a Lorentz metric. Moreover, we consider a torsion-less, symmetric, 
and metric compatible affine connection, namely the Levi-Civita one. The Christoffel 
symbols write as 

- dxg^y) ( 2 . 1 ) 

while the Ricci scalar becomes 

R = g^^idxV^^, - ( 2 - 2 ) 

The four dimensional action that describes F[R) theories then writes as 

5 = ^ / d^x^F{R) + Sm, (2.3) 

with = SttG is the gravitational constant, and where we have also considered the action 
of the matter sector Sm- In the context of the metric formalism, by varying action (2.3) 
with respect to the metric tensor g^^ we obtain the equations of motion: 

F'{R)R^,{g) - \F{R)g^, - V^V,F'(R) + g^,uUF'{R) = (2.4) 

In the above expression and throughout the paper, the prime is assumed to denote differ¬ 
entiation with respect to the corresponding argument (hence in the present case we have 
F'{R) = dF[R)/dR). In addition, stands for the matter energy-momentum tensor 
arising from Sm- 

In order to investigate the cosmological implications of the above gravitational theory, 
we consider as usual a flat Friedmann-Lemaitre-Robertson-Walker (FRW) metric, with line 
element 

ds^ = —dt^-I-a^(t) ^ dx^. (2.5) 
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In this case, the Ricci scalar is calculated to be 


R = Q{2H'^ + H), 


( 2 . 6 ) 


where H = a/a stands for the Hubble parameter and dots indicate differentiation with 
respect to the cosmic time t. Hence, the field equations (2.4), in the case of FRW geometry 
give rise to the Friedmann equation, namely 


-18 


AH{tfH{t) + H{t)H{t)] F"(R)+3 \H\t) + H{t)] F'{R) - ptot = 0, (2.7) 


with ptot the total energy density of the matter fields, and where the Ricci scalar is given 
in (2.6). 


2.2 Superbounce Reconstruction from F{R) Gravity 

The superbounce cosmological scenario was studied in [11], both in a supergravity and non- 
supersymmetric framework, with the most appealing attribute of this scenario being the 
ekpyrotic contraction phase. The ekpyrotic contraction was firstly introduced and studied 
in [12], along with a thorough study on cosmological perturbations. The purpose of this 
Section is to provide an F{R) description of the superbounce cosmological scenario, and 
in order to achieve this we will make use of the reconstruction technique firstly developed 
in [65]. 

The superbounce scale factor is given by [11] 

a(t) ~ (-t + (2.8) 


with t* being the big crunch time, and c a parameter constrained to c > \/6 [11]. Thus, 
the corresponding Hubble parameter is 

= THihry p'**' 


The reconstruction technique developed in [65] is an exact method based on the introduc¬ 
tion of a new variable in place of the cosmological time t, namely the e-folding number N, 
related to the scale factor a{t) through 


e 


-N 


Qq 

a 


( 2 . 10 ) 


Using the new variable N, we can re-express the first FRW equation (2.7) in terms of the 
e-fold parameter N as 


-18 [4H^{N)H'{N) + H^{N){H'f + H^{N)H''{N)] F"{R) 

+3 [H^{N) + H{N)H'{N)] F'(R) - + K^p = 0. (2.11) 

We mention that in relation (2.11), and in the rest of this subsection, the Hubble parameter 
is assumed to be a function of the e-folds N and in addition any derivative appearing in 
(2.11) is defined with respect to N. 
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We now introduce the function G{N) = and therefore we have that 


R = 3G'{N) + 12G{N). 


( 2 . 12 ) 


This relation is very convenient, since it allows us to determine N(R). In particular, the 
Hubble parameter (2.9) can be written in terms of the scale factor as 



and thus eliminating a in favor of N using (2.10) we obtain 


G{N) = 


(2.13) 


(2.14) 


where we have set A = ■ Hence, by combining relations (2.12) and (2.14), we can 

acquire the e-fold parameter as a function of R as 



R 

3A{4 - c2) 


Inserting these in the Friedmann equation (2.11), we can re-express it as 


(2.15) 


-9G{N{R)) [4G'{N{R))+G”{N{R))] F"{R) 




3G{N) + -G'{N{R)) 


F'{R)-^^ + K^ptot = 0, (2.16) 


with G'{N) = dG{N)/dN and G”{N) = d^G(A^)/dA^^. The last step is to express the total 
matter energy density ptot in terms of N{R). Assuming the various matter components pi 
to be independently conserved, i.e. pi + 3F[{1 + Wi)pi = 0, with Wi their corresponding 
equation-of-state parameters, we obtain 

Ptot = (2.17) 


Hence, inserting all the above in the Friedmann equation (2.16), we obtain a differential 
equation for F{R), namely 


aii7^ 


d^F{R) 

di?2 


-|- a2R 


dF{R) 

dR 


F{R) 

2 


F^^SiR 


3(l+^i) 


= 0 , 


where we have set 


and 


ai 


02 


4 — 

2-c 2 
2(4-c2)’ 


Si = 


^ PiO^Q 


[3A(4-c2)]- 


3(l+^-i) 


(2.18) 


(2.19) 


( 2 . 20 ) 
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The solution of (2.18) provides the exact F{R) form that produces the superbounce evo¬ 
lution (2.8). 

Let us first consider pure F{R) gravity, with no matter fluids present. In this case 
the differential equation (2.18) becomes a homogeneous Euler second-order differential 
equation, with solution 

F{R) = ciRP^+C2R^^, (2.21) 

where ci,C 2 are arbitrary parameters, and pi and p 2 are equal to 


Pi 

P2 


-{a2 - oi) -I- \/(o2 - ai)2 -I- 2ai 
2ai 

-(02 - ai) - V(a2 - Qi)^ + 2ai 
2ai 


( 2 . 22 ) 


Finally, it is worthy to examine if this F{R) gravity can become Einstein gravity plus 
curvature terms under some parameter values. This indeed occurs when c S> 1, in which 
case the F{R) gravity takes the form, 


F{R) ~ ciR + C2R~^^^. 


(2.23) 


We now come to the more physically interesting case where matter helds are present. 
In this case, the general solution of (2.18) is easily found to be 


F{R) = 


C 2 P 1 


cipi 




P2 P2{P2 - Pi + 1) 

Y BiC2R^^+P^ + ciRP^ + C2RP^ 


Cl Si 


_P2{^i + 2 -|- P2 ~ Pi) 


P2 


(2.24) 


where ci, C 2 arbitrary parameters, and where we have set 

+ Wi 
C2 

Bi = (2.25) 

P2^i 

It is interesting to note that for 1 the above F{R) gravity, responsible for the super¬ 
bounce generalization in the present pf matter fields, becomes 



- 2c2 


— /52 + 2 


F{R) c^R + aR^ + + A, 


(2.26) 


where we have set C 2 = 1, a = y — 2 -|- 2ci ^ and A = — Aj. 

The F{R) gravity (2.26) can describe both early and late time acceleration, since at 
early times is approximately equal to i? -|- aR^ (Starobinsky inflation) and at late times 
R — \ Por a detailed study of such F{R) forms, and the phenomenological 

determination of the involved parameters, see [70, 118]). 

Let us close this subsection by a remark. In reference [114] a different reconstruction 
procedure was used in order to obtain, in the large and small curvature limits, the F{R) 
gravity that produces the superbounce cosmological solution. The resulting F(R) gravity 
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in the large R limit is quite similar to the one we obtained here, namely of the form R+aR^. 
This was expected to some extent, since the two reconstruction methods overlap at certain 
limits. This overlap however, is in some cases, model dependent and it is strongly affected 
by the particular form of the given Hubble parameter, and in particular from the limit of 
the Hubble parameter at large and small cosmological times. Later on we shall discuss on 
this issue in more detail. 


2.3 Loop quantum cosmology ekpyrotic scenario reconstruction from F{R) 
gravity 


One quite appealing alternative to standard inflation is provided by the ekpyrotic scenario 
[18-22], in which scale invariant perturbations are generated before the Big-Bang phase. 
A recent refinement of the ekpyrotic scenario was presented in [44], in which case LQC 
modifications to the original ekpyrotic scenario are taken into account [35-44]. In the 
context of LQC, the ekpyrotic scenario is realized by using the scalar field potential [44] 


V{^) = 


yggvT6iFGp?i 


1 3pVo ^^/IGwGpd) 

^ ^ 4p,(l-3p)® 


1 2 ’ 


(2.27) 


with the parameter p taking values 0 < p <C 1. As we shall explicitly show in the following 
Sections, the fact that p <C 1 has particularly appealing consequences with respect to early 
and late time cosmological phenomenology, that are absent in other bouncing solutions, 
as for instance in the matter bounce scenario. Since this issue is very important, we shall 
thoroughly discuss it in the end of this Section. 

Before proceeding to the quantitative investigation of the reconstruction procedure, 
we need to make some important clarifications. The potential (2.27) is crucial for our 
analysis, since it guarantees the existence of the ekpyrotic contraction itself. However, 
this potential was derived under the assumption that a phase of ekpyrosis actually exists 
in loop quantum cosmology. Nevertheless, as analytical and numerical investigations in 
the context of loop quantum cosmology indicate [119, 120], for a Steinhardt-Turok type 
ekpyrotic potential, a phase with equation-of-state parameter re S> 1 may not even exist, 
and consequently there may be no ekpyrosis. The reason for this is traced to the behavior 
of the moduli field’s kinetic energy. Particularly, in loop quantum cosmology during the 
contracting phase, and as the moduli field’s kinetic energy approaches Planckian energies 
on the way towards the classical singularity, the Universe bounces. Hence, the Universe 
may never get a chance to reach the rc S> 1 regime, and thus it will bounce in the regime 
where re ~ 1 [119, 120]. Therefore, it is possible that for a large range of parameters in 
loop quantum cosmology, for the Ekpyrotic potential of Steinhardt-Turok type, it is rather 
difficult to obtain ekpyrosis. For an extensive and informative discussion and analysis of 
these issues the reader is referred to Refs. [119, 120]. 

Having these in mind, in the following we assume that an ekpyrotic phase, described 
by the potential (2.27), actually exists and thus in the LQC ekpyrotic scenario the corre¬ 
sponding scale factor and the Hubble parameter are equal to 

a{t) = {aof + , H{t) = (2.28) 
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with Pc the critical density and oq = It is the purpose of this subsection to find 

the F{R) gravity that generates the cosmology described by relations (2.28). As we shall 
demonstrate, we can obtain the cosmology (2.28) without the need for a potential, hence 
modified gravity provides the theoretical framework which can describe ekpyrosis in the 
situation when the corresponding potential is questionable, or does not work completely 
successfully. 

Following closely the lines of the previous subsection, in the case at hand the Hubble 
parameter as a function of the scale factor is given by 

= 4p^ao {a'p — ap'^ , (2.29) 


and by using (2.10) and also the function G{N) = we obtain 

'4 


G{N) = A 


6JV 4JV 

-pcC p — e p 


(2.30) 


where for simplicity we have set the initial scale factor equal to one and also A = 4p^ao. 
Thus, using relations (2.12) and (2.30) we can express the e-fold number as a function 
of the Ricci scalar R as 


iV = ^ln||6Af(3 + 2p) [-2,2A{1 + pf - 9Rp{2> + 2pf + 


-1 


XI - 822/3^2(1 + pf + 4Ai(1 + p) [-32A(1 + pf - 9Rp{3 + 2pf + Mi] 


_2V3hI [-32A(1 + pf - 9Rp{3 + 2pf + ||, 

with Ml = 3y^ Rp{3 + 2pf [64A(1 -|- pf -|- 9Rp{3 -|- 2pf]. 

By substituting expression (2.31) into equation (2.11), we acquire 

, , c^.d?F{R) , , „,dF(R) F{R) 

^l(P) R )—-^ ^2(P) R) —77^-7-f f'm — 0 


(2.31) 


di? 


(2.32) 


with Pm encompassing the contribution of all matter fluids to matter-energy density, and 
where the coefficients Ai{p,R), i = 1,2 read 

Ai{p,R) = - \a32A^p^{ 3 + 2pfMf^^] ^ ^ p)M^^^ + 2^/3jp/2/3l^ 


X |l28A^(l -F pf - 12A(1 + p)A^Mi - 3A^Mf2M2f/^ 
+A^{2M2f/^ [9Rp(3 + 2pf - 8(1 + pf{2M2f/^ 

+4A3(1 + p) ^9Rp{3 + 2pf + 8(1 + pf{2M2f^^^ } 

x|64A^(l + p)3(7 + 4p) -6A(7 + 4 p)A2Mi -3A2Mi(2M2)^/3 

+A^{2M2f^^ [9Rp(3 + 2pf - 8(1 + pf{2M2f^^ 

+2A3 [9Rp(3 + 2pf{7 + Ap) + 16(1 + p)3(2M2)i/3j ^ | 


(2.33) 
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and 


M{R,p) 


72A^p{3 + 


2^^I^A\l + pf 


AA{1 + p)My^ + 


x^2A^[l2 +p{l‘i + 4.p)][Z2A{l+pf + mp{‘i + 2pf -ZMi] 

-3(3 + p)A^Mi{2M2f/^ + 32^3(1 + pf{2, + p)(2M2)^/^ 

+A^{3 + p){2M2)^/^[9Rp{3 + 2pf - 8(1 + p)\2M2)^/^]^ (2.34) 


with Ml = 3^jRp{3 + 2pY [64^(1 + p)'^ + 9Rp{3 + 2pY] and M 2 = 3A^Mi — 32A^(1 + 
pf -9A^Rp{3 + 2pf. 

Solving explicitly the differential equation (2.32) is a challenging and rather formidable 
task, and thus in order to proceed we shall examine certain limits. Particularly, we shall 
study the large and small cosmological time limits of the scale factor (2.28), and by apply¬ 
ing the same reconstruction technique we shall investigate which F[R) gravity generates 
corresponding limits. Note that the large and small time limits correspond to small and 
large curvature limits respectively. 


2.3.1 Large-time approximation of LQC ekpyrotic scenario 

Let us start with the large-time limit. In this case the scale factor and the Hubble parameter 
of (2.28) become 

a{t) K. aoF, H{t) « (2.35) 

and therefore the function G{N) = writes as 

n 2JV 

G{N) = p^e~M. (2.36) 


Inserting (2.36) into equation (2.12) and solving with respect to R we obtain 



6p(2/9 - 1) 
R 


(2.37) 


Hence, substituting (2.37) into equation (2.16), we finally acquire the following differential 
equation: 

( 2 ^) + E = 0’ (2.38) 


with 


Si = [Qaop^{2p - 1 )] 

_ 3(1 -F Wi)p 


Z{l+w^)p 

2 


(2.39) 


The differential equation (2.38) is the Euler non-homogeneous differential equation 
(2.18) which we came across in the previous subsection, thus in pure F{R) gravity (disre¬ 
garding all matter fluids), we get the following solution: 


F{R) = ciRP^ +C2RP\ 


(2.40) 
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with Cl, C 2 arbitrary parameters and with pi and p 2 given by 


3 — P — \/l + lOp + 3 — p + + lOp + (-r,/I1^ 

Pi =-4, P2 = -1-• (2.41) 

Since by construction in the context of the LQC ekpyrotic scenario it is required 0 < p <C 1, 
the solution (2.40) is approximated to 

F{R) = ciR + C2R^/^. (2.42) 


By taking into account the presence of ordinary matter fields, and using the technique 
we employed in the previous subsection, the reconstructed F{R) gravity which is solution 
of (2.38) reads 


F{R) = Cl 

i 

-E 


Si 


P2{oii — Pl) 
SiC2 
P2{oii — P 2 ) 


- 


ClPl 


P2{P2 - Pl + l) 


+ ciRP^ 


P2 


(2.43) 


and therefore for 0 < p <C 1 it becomes 

F{R) = ciR + ^( ~ ) R^ + C2R^/‘^ + A, (2.44) 

with A = X^j 5 'i(c 2 — 2ci). Thus, by looking relations (2.44) and (2.42), we deduce that 
the effect of the matter fluids is to introduce a cosmological constant plus subdominant 
curvature corrections. This result is quite similar to the reconstructed F{R) which pro¬ 
duces the superbounce solution given in relation (2.26), but now the late-time behavior 
is different. However, note that relation (2.26) gives the superbounce generating F[R) 
gravity for all values of the Ricci curvature R, whilst (2.44) is valid only for small values 
of the curvature (large cosmological times). Thereby, the result (2.26) has a larger validity 
region and moreover is phenomenologically more appealing. 


2.3.2 Small-time approximation of LQC ekpyrotic scenario 

In the small-time region the LQC ekpyrotic scale factor and Hubble parameter (2.28) 
approximately becomes 

a{t) « 1 -h H{t) K. aopt (2.45) 

and therefore, using relations (2.10) and (2.12), the e-fold parameter N can be expressed 
as a function of the Ricci scalar R through 


N = \n 


( 24ao/0 -l- R 
V 30ao/0 


(2.46) 


Therefore, the differential equation (2.11), without taking into account any matter fluids, 
can be cast in the following way 


2 (72ao/o^ - I 800 / 9 R - R^) F"{R) + 3(4aop + R)F'{R) - 5F{R) = 0. (2.47) 
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This differential equation can be solved using hypergeometric functions, namely 


F{z) = 2 ^ 1 ( 0 ;,/3, 7 ; z) 


(2.48) 


where 


z = 


6ao/9 — R 
30aop 


with Ai = Gaop and A 2 = — 24ao/o, and where a, /3 and 7 satisfy 

5 3 1 

al3 = -, a + P + 1 = 7 = - 2 j 


which accepts two sets of solutions, namely 


1 


a = - -5 ± V65 , /3 = - -5 =F V65 


(2.49) 


(2.50) 


(2.51) 


Since we are in the large curvature limit, we can expand the above general solution in terms 
of l/R, obtaining 

+3600/9^ {a(l + a)r[—a + /3]r[7]} 

[(1 — /3 + 7)(7 - /3) - 8(1 + a - 7)(7 - /3) + 16(1 + a - 7)(2 + a - 7 )] 


2{1 + a - (3){2 + a - f3) r[/3] r [7 - a]R^ 




(2.52) 


Hence, given the values of a and /3 from (2.51), we deduce that the most dominant term 
is the first one, namely 




(2.53) 


Let us make some comments here. In [114] it was shown that the large-i? reconstructed 
gravity which generates the LQC ekpyrotic scenario is of a different form compared to 
the above relation. The source of the difference is that in the present work we use a 
more accurate limiting behavior for the small-time scale factor (namely relation (2.45)) 
comparing to [114]. In particular, apart from the method we used above, there is another 
well-known reconstruction, firstly developed in [ 66 ], where the reconstruction is achieved 
using an auxiliary scalar field. In reference [114] the auxiliary-field reconstruction was 
used in order to face the same issue, namely to extract F{K) gravities that produce the 
superbounce and LQC ekpyrotic scenario. Concerning the superbounce with ordinary 
matter fluids present, both in the large and small curvature limits, the results are the 
same to the ones obtained in this work. On the other hand, concerning the LQC ekpyrotic 
expansion, in the small-curvature limit the two methods give the same results, while in the 
large-curvature limits a difference appears. The reason behind this difference is traced in 
the particular form of the scale factor and Hubble parameter and also their limits when t 
tends to zero, a case which corresponds to large curvatures. Specifically, the scale factor of 
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the LQC ekpyrotic scenario given in relation (2.28), when t tends to zero is actually equal 
to one, i.e. limi^o®(^) = Ij and moreover the limit of the Hubble parameter (2.28) is also 
zero, namely lim^^o H(t) = 0. This is the result that is taken into account in the auxiliary 
field reconstruction technique and is materialized explicitly in reference [114], resulting to 
a different F{R) gravity in the large curvature limit (see also [121] for similar results). 

On the other hand, in the approach we adopted in this article, if we consider the 
above limits for the scale factor and Hubble parameter, then the reconstruction technique 
yields trivial results or simply does not work, because G{N) would be identically equal to 
zero. Thus, we need an exact functional dependence of the scale factor and of the Hubble 
parameter, with respect to cosmic time t, in order for the method to work perfectly. That 
is why we expanded the scale factor in power series of t and we kept the leading order 
terms as in (2.45). In effect, we could say that in this particular case the auxiliary held 
method is a bit more accurate, to leading order at least. 

In summary, the above artihcial difference is traced on the particular form of the 
scale factor and of the Hubble parameter, and if someone was able to solve the differential 
equation without taking any limiting cases (a rather formidable task), then there would be 
no difference between the two approaches. An exemplihcation to support this argument is 
the superbounce case which we were able to solve analytically, in which case we had absolute 
concordance. Hence, in principle, both reconstruction techniques lead to the same results, 
but extra caution is required if the differential equation that yields the reconstructed F{R) 
gravity, cannot be solved analytically. 

Before closing this subsection an important remark is in order. As we demonstrated 
above, the loop quantum ekpyrotic cosmology can be described by pure metric F{R) grav¬ 
ity, with pure referring to the absence of matter. There exists a very well developed strong 
argument, according to which the full loop quantum evolution is successfully described by 
Palatini F{R) gravity [122] and not by metric F{R) theories. Notice however that a metric 
F{R) gravity is mathematically equivalent to Palatini F{R) theory in the absence of mat¬ 
ter fluids, i.e. for pure vacuum solutions, as those we presented in the previous analysis. 
It would be interesting to study the description of the loop quantum ekpyrosis, when it 
exists, in the presence of matter helds, by using both metric and Palatini F{R) gravity. 
We hope to address this issue in a future project. 

3 Superbounce and loop quantnm ekpyrotic cosmology from F{G) grav¬ 
ity 

In this Section we will reconstruct the F(G) form that can lead to the superbounce and loop 
quantum ekpyrotic cosmology realizations. Before analyzing the reconstruction procedure, 
in the following subsection we briefly review the basic features of F{G) gravity and we 
describe the geometric background we will use in its cosmological application. For a detailed 
account on these issues, see [54] and references therein. 
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3.1 F{G) gravity and cosmology 

Let us review briefly the F{G) gravitational modification and its cosmological applica¬ 
tion. This class of modified gravity is based on the use of the Gauss-Bonnet combination 
G = + R^ypfjR^'^P^, with Rp,upa and Rpy the Riemann and Ricci tensors 

respectively [74-78]. In particular, one can construct gravitational modifications using an 
arbitrary function F{G), which prove to lead to interesting cosmological behavior [79-87]. 
The F{G) action takes the form [79, 80] 

5 = ^ / [R + F{G)] + Sm, (3.1) 

with = 8itG is the gravitational constant, and where we have also considered the action 
of the matter sector Sm- The corresponding field equations for a general metric write as 

Rf.. - lgt.uF{G) + (2RRp, - ARppRP + 2RP^^R,p^r - R^aupRpa) F'{G) 

+4 [VpV.F'iG)] RPp - igp, [VpV^F\G)] RP^ + 4 [VpV<,F'(G)] g^Pgf^^Rpau0 
-2 [V^V,F'(G)] R + 2gp, [□F'(G)] R 

-4 [□F'(G)] Rp, + 4 [V^V,F'(G)] RP = (3.2) 

where is the matter energy-momentum tensor arising from Sm- In the case of the FRW 
metric (2.5), the above equations give rise to the two Friedmann equations 

+ F{G) - GF'{G) + 2dH^GF''{G) = 2K^pm (3.3) 


AH + + F{G) - GF'{G) + 16HG (^H + F"(G) 

+8H^GF''{G) + 8H^G^F'''{G) = -2K^pm, (3.4) 


where pm and pm are the matter energy density and pressure respectively. Note that in 
FRW geometry, the Ricci scalar is given by (2.6), namely R = 6{2H^ + H), while G reads 

G = 2AH^ (^H + . (3.5) 


3.2 Ekpyrotic scenario reconstruction from F{G) gravity 

We now use the reconstruction method developed in Refs. [84, 85], for finding the F{G) 
gravity, in the large cosmic time limit, which corresponds to the late-time era of the ekpy¬ 
rotic scenario. For an alternative method to the one we shall employ, see Refs. [86, 87]. 
We start from the pure F{G) gravity action. 


s = I ^ + 

and by introducing the auxiliary scalar field (p the above action becomes 


S 



V{(t>)-mG 


(3.6) 


(3.7) 
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Varying the action (3.7) we obtain 


V'(<^)+e'(</>)G = 0, 


(3.8) 


and if this equation is solved with respect to (p we obtain the functional dependence of 
the auxiliary scalar as a function of G, namely (j) = </>(G). Then, by substituting into the 
action (3.7) we obtain the function F{G), which is 

F{G) = -V{cj,{G))-mG))G. (3.9) 


Hence practically, if one finds 4>{G), the F{G) gravity follows straightforwardly. For the 
flat FRW background (2.5) we obtain the following field equations; 

= 0, (3.10) 

dt 

^ (2H + 3H‘^] - VU) 16HH^ - 16H^^ = 0, (3.11) 

\ J dF dt dt 

and upon combining them we obtain 

+ = 0. (3.12) 

dF dt dt 

Solving Eq. (3.12) with respect to ^(t) we acquire 


N 1 /"* , a(x) . 

(3.13) 

with IF (t) being equal to 

W{t) = 2 / dx 

J a{x) 

(3.14) 

and with the prime denoting differentiation with respect to x. Combining (3.12) and (3.13) 
we obtain 

V(t) = ^Hitf - 3a(t)H(t)W(t). (3.15) 

Having the functions F(t) and ^(t) at hand we can easily obtain the F{G) gravity by 
using Eq. (3.9). The large t limit of the scale factor is a{t) = and the corresponding 

Hubble rate is F7(t) = j, therefore the function ^{t) reads as 

~ 8F^p{l + p) ’ 

(3.16) 

and the function V (t) is equal to 


V{t) = F 

Fk^ Fk^p + 1) 

(3.17) 

By substituting into Eq. (3.8), and upon solving with respect to t, 

we obtain 

23 / 431/4 (/_p3)l/4 
^ “ GV4 

(3.18) 
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and finally substituting this into Eq. (3.9), we obtain the final form of the F{G) gravity 
in the large cosmic time limit, namely 


F(G) = - 




ii‘^p{l + p) 

From this relation we can see that the resulting behavior at large t is rsj Vg. 


(3.19) 


3.3 Superbounce Reconstruction from F{G) Gravity 

Let us now reconstruct the F{G) form that leads to the superbounce realization. It proves 
convenient to use the well-known reconstruction technique from [84, 85], which we briefly 
described in the previous section. Following the same steps as in the previous section, we 
may easily obtain the functions ^{t) and V{t) of Eqs. (3.13) and (3.15). Using the Hubble 
rate and the scale factor from Eqs. (2.9) and (2.8), we easily obtain the function V{t) 
which reads 

12 (c^ - 2) 


V{t) = - 


while the function ^{t) reads. 


m = 


(2 -I- c^) (U - 
cH{t — 2U) 


(3.20) 

(3.21) 


32(2 + c2) ■ 

Having these at hand, we easily solve (3.8) with respect to t, and we obtain the following 


solutions: 


2 


t — U 


■(c2-2)k2g' 

1/4 

2 23 / 431/4 


1/4 

c® 

c® 



5 t — “h 




(3.22) 


Hence, by substituting these solutions into (3.9), we easily obtain two solutions for the 
F{G) gravity that generates the superbounce, which read 


Fi{G) = 


,4 c 

^ F2{G) = 




32 (2 c2) 


32 (2 c2) «:2 


(3.23) 

We need to note that the method of Refs. [84, 85] is much more direct in comparison to 
the method of Refs. [86, 87], since in the latter case the procedure is more complicated. 


4 Superbounce and loop quantum ekpyrotic cosmology from F{T) grav¬ 
ity 

In this Section we will repeat the above reconstructions in the case of F{T) gravity, i.e. we 
desire to reconstruct the F{T) forms that can lead to the superbounce and loop quantum 
ekpyrotic cosmology realizations. We start our analysis by a brief introduction to F{T) 
gravity. 
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4.1 F{T) gravity and cosmology 

In general, in the torsional formulations of gravity, one uses as dynamical fields the vierbeins 
(we use Greek indices for the coordinate space and capital Latin indices for the tangent 
spacetime), which at each point of spacetime form an orthonormal basis for the tangent 
space. The metric tensor is obtained from the dual vierbein through 

Qt^uix) = riAB e^{x) ef (x), (4.1) 

with r]AB = diag(l, — 1, — 1, — 1). Contrary to General Relativity and its curvature-based 
extensions {F{R), F{G) etc), where one uses the torsionless Levi-Civita connection (2.1), in 

the present framework we use the curvatureless Weitzenbock connection defined as = 
e\ dfj_e^ [123]. Hence, all the gravitational information is embedded in the torsion tensor, 
written as 

w-^ w-^ 

Tl, = {d^et - d,ei). (4.2) 

Similarly, the contorsion tensor reads as = — i^T^p — T^p — Tp , while we addi¬ 
tionally define -|- 6p T°‘’^ — Sp One can now construct the torsion 

scalar through contractions of the torsion tensor (similarly to the construction of the Ricci 
scalar from the Riemann tensor) as [124-126] 

T = ^r^^'Tpp, + ir^'^'r.pp - (4.3) 

In the simple, teleparallel equivalent of General Relativity, the Lagrangian is just T. 
However, inspired by the F{R) extensions, one can extend T to a function F{T), and thus 
the action in the F{T) modified gravity writes as [88-113] 

R = ^ / d^xe [R(r)] + Sm, (4.4) 

where e = det(ej^) = y/—g, = 8ttG is the gravitational constant, and we use units 

where the light speed is equal to 1. Note that we have also added the matter action Sm- 
Therefore, varying the action (4.4) with respect to the vierbeins leads to the held equations: 

e-'ap(ee^5p'^")Rr + e^Rp'^"5p(r)RTr - FTe\T%xSp^^ + -e\F{T) = -e^A^p^, (4-5) 

m 

with Ft = dF/dT, Ftt = d'^F/dT'^, and Tp'^ denoting the usual energy-momentum tensor. 
Note the great advantage of F(T) gravity, namely that its held equations are of second 
order, while in F{R) gravity they are of fourth order. 

In order to apply F{T) gravity in a cosmological framework, we have to consider an 
ansatz for the vierbeins that corresponds to the FRW metric (2.5), namely 

Cp = diag(l,a,a,a). (4.6) 
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In this case, the field equations (4.5) give the two Friedmann equations as: 


TFt 

3 


„ + „ Pm — 0 
D 3 


H(^Fx + 2FFxt) 


2 {Pm + Pm)j 


with Pm and pm the matter energy density and pressure respectively, and where 


(4.7) 

(4.8) 


T = 


(4.9) 


as it arises from (4.3) for the FRW ansatz (4.6). Hence, the General Relativity equations 
are obtained for F{T) = T. 


4.2 Superbounce Reconstruction from F{T) Gravity 

Let us now suitable reconstruct the F{T) form that can gibe rise to the superbounce 
realization [11], i.e. to the scale factor (2.8), namely 

a(t)-(-t + G)2/^', (4.10) 


with the big crunch time, and c > y/d the ansatz parameter. As usual we use the e- 
folding number N through e~^ = and thus for the Hubble parameter we obtain H = 
2a~'^ /^/2. The great advantage of F{T) reconstruction, comparing to F{R) reconstruction 
of subsection 2.2, is that the torsion scalar T is related straightaway with the Hubble 
parameter through (4.9), while the Ricci scalar is related to both the Hubble parameter 
and its derivative. Hence, in the present case we can immediately express N as a function 
of T as 



(4.11) 


with 


7 24 _,2 

A - ^4 «0 ■ 


(4.12) 


Additionally, using the matter energy density conservation, we can express pm as 


Pm — 


Pi0% 




(4.13) 


Inserting these in the first Friedmann equation (4.7) we obtain 


rj. dF{T) 

dr 


HT) 

2 


+ J2S.T 


3(1+Wi) 

1 - 


= 0 , 


(4.14) 


with 


Si 


_-3(l+«)i) 

^ PiOOiQ 




(4.15) 


The solution of (4.14) provides the exact F{T) form that produces the superbounce evo¬ 
lution. 
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We can immediately see that (4.14) is much more simple than the corresponding equa¬ 
tion (2.18) of the F{R) case. In particular, in the case of absence of matter fields, the 
solution reads 

F{T) = Cl i-T)^ , (4.16) 

where ci is an arbitrary parameter (note that the presence of the minus sign in front of T 
was expected, since T = —6H^ < 0 in the conventions used in this work). Hence, simple 
teleparallel equivalent of General Relativity, without matter fields, cannot give rise to the 
superbounce evolution. However, when matter fields are present, the corresponding F{T) 
is found to be 

1 «—T 3(l+u)i) 

F{T) = ci{-T)-2 , (4.17) 


where 


ZK C PiOCLq 

6(1 -k Wi) - 


(4.18) 


As we observe, for c = 3(l+'u;j) (c = 3 for simple dust matter) we find that the superbounce 
is realized in the case of standard teleparallel equivalent of General Relativity, i.e. standard 
gravity, plus corrections. Finally, note that the above F{T) form is in general different than 
the one required for matter bounce realization in F{T) gravity [127]. 


4.3 Loop quantum cosmology ekpyrotic scenario reconstruction from F(T) 
gravity 

Let us now reconstruct the F{T) form that gives rise to the LQC ekpyrotic scenario (see 
also [128] for a different approach). As we mentioned in subsection 2.3, the desirable scale 
factor and Hubble parameter should have the form (2.28), namely 

a{t) = {aot^ + , Hit) = (4.19) 

with 0 < p <C 1, and where pc is the critical density and oq = The Hubble parameter 

can then be expressed as 

= 4p^ao — a7^ , (4.20) 

and since in H(T) gravity T = (relation (4.9)), we can immediately find that 

T = —24p^ao — a7^ . (4.21) 


Therefore, introducing the e-folding number N = ln(a/ao) and eliminating a we can express 
N as a function of T through 


NiT) = In 



+ 


(f)* 



(4.22) 
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with M 2 = 2 - Finally, inserting N{T) from (4.22) into pm- 

equation (4.13) and then into the Friedmann equation (4.7) we obtain 

3p(l+^i^j) 

2 

= 0, (4.23) 

with 

Zp{1 + Wa) rt 

Xi = 3^^K^p^o. (4.24) 

The solution of the above differential equation will give the F{T) form that generates the 
LQC ekpyrotic scenario. Unfortunately, this equation cannot be solved analytically, and 
hence, similarly to subsection 2.3, we solve it in the large and small time regimes separately. 


,dF(r) F{T) 
“dr ^ 






4.3.1 Large-time approximation of LQC ekpyrotic scenario 

At large times the required scale factor (4.19) becomes a{t) ^ , and thus = 

p^ (oo/a)^^^. Thus, the N{T) expression is significantly simplified, and it reads 


N{T) = ^\n 



(4.25) 


Hence, inserting N{T) from (4.25) into pm-equation (4.13) and then into the Friedmann 
equation (4.7) we obtain 


^ dF(r) 

dr 


m+Y.yd-Tr=o, 

i 


(4.26) 


with 


0/1 I N 3 (l+u;i)p 

hi = K^Pioa^ (6p^) 2 


Oti = 


3(1 + Wi)p 


Equation (4.26) can be easily solved as: 


r(r) = ci(-r)^-^ 


2 Ti 


2 aj - 1 


{-TT , 


(4.27) 


(4.28) 


which is the form of F{T) gravity that can generate a LQC ekpyrotic scenario at large 
times. 


4.3.2 Small-time approximation of LQC ekpyrotic scenario 

At small times the required scale factor (4.19) becomes a{t) ~ 1 + thus H{t) ~ 

ttQpt. Thus, the N{T) expression becomes 


N{T) = In 



12 pal) 


(4.29) 
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Therefore, inserting N{T) from (4.29) into pm-equation (4.13) and then into the Friedmann 
equation (4.7) we obtain 



= 0 


(4.30) 


with 


ry 2 

^ PiO 

ai = -3(1 + Wi). 


(4.31) 


Equation (4.30) has the solution 


F{T) = ci{-T)^ +2Y,Zi2Fi 



(4.32) 


with 2 Fi{a,b-, c; z) the hypergeometric function. Hence, this is the form of F{T) gravity 
that can generate a LQC ekpyrotic scenario at small times. Finally, since at early times 
we are in the large torsion regime (T 3> 1), we can approximate the above solution as 





(-r)i + {-rr > • (4.33) 


FiT)^ci (-T)^+2j^Z, 


5 Stability of F{R), F{G) and F{T) superbounce and LQC ekpyrosis 

In the previous Sections we investigated the superbounce and LQC ekpyrosis in the frame¬ 
work of F{R), F{G) and F{T) modified gravities. A crucial issue in every cosmological 
evolution is whether at the perturbation level it exhibits instabilities, which thus could 
constrain or exclude the scenario. In this Section, we perform such stability investigation, 
by linearly perturbing the above solutions. We mention that we are interested in the pure 
gravity case, in order to obtain a first picture for the behavior. The study of the full 
perturbation behavior, including the matter fields, lies beyond the scope of this work, and 
will be addressed in a separate project. 

5.1 Stability in the F{R) reconstructions 

In the context of the reconstruction method [65] , the stability study was performed in detail 
in [129], the notation and formalism of which we adopt in this Section. The stability of 
the solutions which we obtained in the previous Sections, can be examined by performing 
linear perturbations of the F{R) solution and specihcally of the auxiliary function G{N). 
We consider its perturbation to be of the form [129] 


G{N) = g{N) + 6g{N) 


(5.1) 


and we insert it in the first FRW equation (2.16). Since the background quantity g{N) 
satisfies exactly equation (2.16), we easily acquire the equation for the perturbation 6g{N), 
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namely 


giN) 


(PF{R) 


dR^ 


R=R\ 


5"g{N) 


+ |3g(JV) [4g'(iV) +s"(JV)] 
+ I 129 (JV) [4g'(JV) + /(AT)] 


R=Ri 


+ 


3g{N) - -g'{N) 


d?F{R) 


di ?2 


R=R\ 


5'g{N) 


R=R\ 


+ [-Ag{N)+2g'{N)+g"{N)] 


d^F{R) 


di ?2 


+ 


1 dF{R) 


R=Ri 3 dR 


R=Ri 


6 g{N) = 0, (5.2) 


with Ri = 3g'{N) + 12g{N). From this equation one may directly extract the stability 
conditions for the perturbations of G{N). In particular, they read 


Ji 


6 [dg'{N)+g"{N)]F'"{R) g'{N) 

F"{R) ^ giN) 


(5.3) 


and 

_ 36[dg'{N)+g"{N)]F"'iR) dg'jN) 3g"iN) F'{R) 

' F"{R) giN) giN) giN)F"iR) 


(5.4) 


Let us now examine whether the super bounce and LQC ekpyrotic FiR) forms found in 
Section 2 satisfy the above stability conditions. Considering the superbounce cosmological 
solution (2.21), without the presence of matter field, the stability conditions (5.3) and (5.4) 
respectively become 


Ji = 6 + - 2 c 2 [3^iciQf (pi - 2)(pi - l)pi + 3^^C2gf (p2 - 2)(p2 - 1)P2] 

X [3P^ciQ{^ipi - l)pi +3P^C2Ql^ip2 - 1)P2] > 0, (5.5) 

with Qi = —g4(c2 — 4)e“'^^'^, and 


J 2 = -12 - 6c^ + 30“^ + (^ciQ^^ ^pi + 02 ^ 2 ^ ^^ 2 ) A ^ 

-12c^Q2 [ciQ^^-^(pi - 2)(pi - l)pi + C2Q^2'~\p2 - 2)(-l + ^ 2 )^ 2 ! Qg ' > 0 (5-6) 


with Q 2 = 3Ae~'^^’^i4: - c^) and Q 3 = ciQ^^~‘^ipi - l)pi + C 2 Q 2 ^~‘^iP 2 - 1)P2- The most 
interesting subcase that we found in Section 2 was when the parameter c is large, in which 
case the above stability conditions become 


= 4c2 + 6 > 0 (5.7) 

J 2 = 30^^ + ISc^ - 12 > 0. (5.8) 


Hence, in the interesting case of the large-c regime, the FiR) superbounce is free of insta¬ 
bilities. 

Considering now the superbounce cosmological solution (2.21), with the presence of the 
matter field, i.e. the FiR) form (2.24), the stability conditions (5.3) and (5.4) respectively 
become 
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= 6 + - 6c2Qi [3^l-"ci(pi - 2)(pi - l)piQ{^-^ + r^-^C2{p2 - 2){p2 - l)/92Qf " 

+ 3^^-2gi(p2 - 1 )P 2 (P 2 + l)Qf" - 3'5^+^2-3g5.+P2-3^^^^^ ^ _ 2 )^ 5 . + 

_j_ 35j+P 2 Ig5,+P2 lg2((5j + P2)((5i +/52 + l)(5i + P 2 + 2 ) 

X 3^1"^Ci(/9i - l)/9i(5f + 3^2"^C2(/92 “ l)/02Qf + 3'’2-lgip2(p2 + l)Qf 

- + P2 - 1)(<5* + /92)Qi^+'’^"' + + P2 + 1)(<5* + P2 + 

(5.9) 

with Qi = —A{c^ — 4)e“'^^^, and 


-1 


J 2 = -12 - Gc^ + 30^^ + ^ I 


,P2-1 


(Qs + <54) ^ + C2P2Q2 

+Qf gi(l + p2) - + P2) + QI^^^^+\2{2 + 5^ + P2)\ } 


+36Ac2e-"'^(c2 - 4) (Q 3 + Qi)~^ 


CiQ ^2 ^(pi - 2)(pi - l)pi + C2Q2^ ^(P2 - 2 )(p2 - 1)P2 

+Q?~^gi(/02 - 1)P2(/02 + 1) - Q 2 + P2 - 2){5i + P2- l)((5i + P 2 ) 


+Q 2 ~^^^ ^(l 2 {^i + P 2 ){^i + P 2 + l)(<5i + P 2 + 2 ) 


(5.10) 


with (^2 = 3^e '^^^(4 - c^), Qa = ciQ^^ ^(/?i - l)pi + C 2 Q 2 ^ ^(p 2 - 1 )P 2 and Qa = 

^2^ 91^2(1 + P2) — Q2^^^ <? 2 ( 5 i + P2 — l)(( 5 i + P2) + Q2^^^ + < 5 i + P 2){2 + 5 i + ^2)- 

Thus, in the large-c regime they become 


Ji « Gc^ (5.11) 

J2-c\ (5.12) 

We mention that by large c we practically mean c > 10 and we also recall that c > y/6 in 
order a superbounce exists. From relations (5.11) and (5.12) it is obvious that Ji > 0 and 
also J 2 > 0, and therefore the superbounce generating F{R) gravity, with or without the 
presence of matter, is stable under the perturbation (5.1). 

Considering the LQC ekpyrotic F{R) gravity that we reconstructed in subsection 2.3, 
we mention that we extracted only approximate solutions in the large and small time 
regimes. Although the consistent perturbation analysis should be performed in the full 
expressions, for completeness we examine the stability of the approximate F{R) form given 
in (2.53). In this case the stability conditions (5.5) and (5.G) become 


and 


■^1 - 1 - 2 « + > 0 


(5.13) 


J 2 — —12 + 


9e^ G0e^(2 + a) 3e^ + 12 (e^ - l) 

5e^ — 4 [e^ — 1) (1 + a) 


e ^-1 


< 0 . 


(5.14) 
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As we can see, for small e-folding values N these conditions are not satisfied and thus the 
approximate solution exhibits instabilities. However, the numerical investigation of the 
exact solution is necessary before we conclude on the stability of the F{R) LQC ekpyrotic 
scenario. 

5.2 Stability in the F{G) reconstructions 

We now examine the stability in the F{G) reconstructions of Section 3, repeating the steps 
of the previous subsection. The general F{G) stability analysis was performed in [85], 
and thus we adopt the notation of that work. Inserting the perturbation (5.1) into the 
Friedmann equation (3.3) in the absence of matter fields, we extract the following stability 
conditions: 

^>0, ^>0, (5.15) 

Ji Ji 

where Ji stands for 

Ji = 288g{NfF"{G) (5.16) 

J 2 = 432g(iV)2 {{2giN)+g'iN))F'\G) 

+8g{N) [g'{Nf + g{N){4g'{N) + g''{N))] F”{G)} (5.17) 

Js = 6 {l + 24g{N) {-8g{Nf + 3g'{Nf + 6g{N)[3g\N) + /(iV)]} F"(G) 

+24giN)[4giN)+g'iN)] {g\Nf + g{N)[4giN)+g^'iN)]}F"{G)]. (5.18) 

In the case of the superbounce F{G) gravity, we found two solutions Fi{G), i = 1,2, 
given in relation (??). Thus, for the function Fi(G) we calculate: 

^ = hc‘^ -2) - l) > 0 (5.19) 

and 

J 3 _ (-f c® -k 108c^ -k 26c2 - 8 ) ^ (lie® - 46c^ -k 30c^ - 4) 

-'^(‘1 - 

with A = ■ Hence the Fi{G) solution is conditionally stable (note however that in 

the large-c region the first condition is not satisfied and thus instabilities appear). 
Similarly, for the case F 2 {G) the stability conditions (5.15) write as 

^ = I (c^ - 2) - 1 ) > 0 (5.21) 

Ji 2 V / 
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and 


J 3 _ 4A^ (16 - 60c^ - igOc'^ + 207c®) 

’^1 Qa 

^ 4 ^ 3 g-c 2 Ar ^_32 200c2 - 152c^ + 122c® - 68 c® + llc^®) 

^ (16 - 92c2 + 30c^ - 46c® + 11c®) ^ (-8 + 26c2 + llTc"^ 

16V6^2^2 + c2)(4 - 4c 2 + c^)v^(2 - Ilc2)(c2 - 2) ^ 


99c®) 


(5.22) 


where Qa = (2 — llc^) 4A^{c^ — 2) + . Hence the i^i(G) solution is conditionally 

stable (note however that in the large-c region the first condition is not satisfied and thus 
instabilities appear). 


5.3 Stability in the F{T) reconstructions 

Let us now discuss on the stability of the obtained superbounce and LQC ekpyrosis in 
F{T) gravity. As we have mentioned, a crucial difference, and a main advantage, of F{T) 
gravity, comparing to F{R), F{G) and other curvature gravitational modifications, is that 
its equations of motion are of second and not of fourth order. This can be immediately 
seen in the Friedmann equation (4.7) of F{T) gravity, in which only H and not H ox H 
appears, and compare it with the Friedmann equations of F[R) and F{G) cases, equation 
(2.7) and (3.3) respectively, where all H, H and H are present. Hence, in order to study the 
stability under linear perturbations in F{T) gravity solutions, we cannot apply the method 
we used in the F{R) and F{G) cases, which was based on perturbing the auxiliary function 
G(A^) = H‘^{N) as in (5.1), since in this case the results will be trivial. The method used in 
the F[R) and F{G) cases is strictly speaking a method of perturbations of the dynamical 
system formed by the FRW equations. In the F{T) case, the stability examination requires 
to perform the complete perturbation analysis, starting from perturbing the vierbeins. This 
analysis has been performed in detail in [89, 90, 104], and therefore we do not repeat it here. 
However, using those results, we can easily deduce that the simple power-law F{T) forms 
of the superbounce ((4.16) and (4.17)) and LQC ekpyrotic ((4.28) and (4.33)) solutions 
are stable (power-law solutions and especially the square root ansatz are very common in 
F{T) gravity [88-113]). 


6 Conclusions 

In this work we investigated the superbounce and the loop quantum cosmological ekpyrosis, 
in the framework of various modified gravities. Bouncing solutions can be an alternative 
to inflation and thus it is both necessary and interesting to see whether it can be obtained 
naturally, for suitably chosen classes of gravitational modification. We focused our investi¬ 
gation in the F{R), F{G) and F(T) modified gravities and in each case we reconstructed 
the forms of modification that give rise to superbounce and LQC ekpyrotic evolutions. The 
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reconstruction methods that were used in this work are not the only ones that one could 
follow. Indeed, depending on the modified gravity that is used, one could use different 
techniques, which usually lead to the same exact results, or similar results when approx¬ 
imations are used. We will perform a careful comparison and discussion of the various 
reconstruction methods in a separate publication. 

Concerning F{R) gravity, we found that the superbounce evolution at the large- 
curvature regime is obtained from R+aR'^, while at small curvatures from R+ciR~^^‘^ +A, 
which agrees with the results obtained in [114] from a different approach. Similarly, the 
LQC ekpyrosis is realized from hyper geometric functions that can be approximated by 
power-law F{R) forms. Concerning F{G) gravity, we found that the superbounce and 
LQC ekpyrosis are obtained from F{G) forms having G, y/G and IjVG factors. Further¬ 
more, in the case of F{T) gravity, we found that the superbounce is generated by power-law 
F{T) forms, while the LQC ekpyrosis from hyper geometric F[T) functions that can be ap¬ 
proximated to power laws. Moreover, performing a linear perturbation analysis, we showed 
that the obtained solutions are conditionally or fully stable. 

It is expected that bouncing cosmology should occur due to quantum gravity effects. In 
this respect, the important step in its realization is the construction of a realistic bouncing 
cosmology within a specific effective gravity theory. It seems that very similar classes 
of different effective gravities, i.e. power-law, polynomials and hypergeometric functions, 
naturally lead to realization of the proposed bounce universe. This indicates that indeed 
quantum gravity, which may give the origin to specific effective F[R), F[G) or F{T) gravity 
theories, may be responsible for the occurrence of bouncing universe. 

Finally, we mention that a crucial issue is to analyze in detail the processing of per¬ 
turbations through the bouncing phase in these scenarios, which would lead to extract 
observable signatures than can be confronted with observations and in particular with the 
Planck data. This is the main test that would lead to constraining or excluding the scenar¬ 
ios at hand, similarly to the various inflationary models. In addition, a very compelling task 
related to the perturbations analysis is the deviation from scale-invariance in the power 
spectrum, which would result to a blue or red tilt in the spectrum of tensor perturbations. 
These necessary analyses he beyond the scope of the present work, and are left for a future 
investigation. 
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